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Abstract — This paper focuses on controlling the absorbing set 
spectrum for a class of regular LDPC codes known as separable, 
circulant-based (SCB) codes. For a specified circulant matrix, 
SCB codes all share a common mother matrix, examples of which 
are array-based LDPC codes and many common quasi-cyclic 
codes. SCB codes retain the standard properties of quasi-cyclic 
LDPC codes such as girth, code structure, and compatibility with 
efficient decoder implementations. In this paper, we define a cycle 
consistency matrix (CCM) for each absorbing set of interest in 
an SCB LDPC code. For an absorbing set to be present in an 
SCB LDPC code, the associated CCM must not be full column- 
rank. Our approach selects rows and columns from the SCB 
mother matrix to systematically eliminate dominant absorbing 
sets by forcing the associated CCMs to be full column-rank. We 
use the CCM approach to select rows from the SCB mother 
matrix to design SCB codes of column weight 5 that avoid all 
low-weight absorbing sets (4,8), (5,9), and (6,8). Simulation 
results demonstrate that the newly designed code has a steeper 
error-floor slope and provides at least one order of magnitude 
of improvement in the low error rate region as compared to an 
elementary array-based code. 

I. Introduction 

It is well known that finite-length low-density parity-check 
(LDPC) codes suffer performance degradation in the high 
signal-to-noise ratio (SNR)/low frame error rate (FER) region. 
This degradation is commonly referred to as the error floor. 
Prior work indicated that certain sub-graphs called trapping 
sets H], and, in particular, a subset of trapping sets called 
absorbing sets ||2) are a primary cause of the error floor. 
Absorbing set is a particular type of a trapping set that is stable 
under bit-flipping decoding. This paper improves performance 
by controlling absorbing sets for a class of regular LDPC 
codes, known as separable circulant-based (SCB) codes, that 
are constructed as an arrangement of circulant matrices. 

Recent papers have proposed methods to lower the error 
floor by improving the absorbing set (or trapping set) spec- 
trum. For example, small trapping sets can be avoided by 
introducing additional check nodes Q, or by increasing the 
girth [4|. The algorithm in [5 ] constructs quasi-cyclic codes 
from Latin squares so that the Tanner graph of the code does 
not contain certain trapping sets. Recent results 151, IfTl, ll8ll 
have also investigated the errors floor of certain practical codes 
in terms of their trapping/absorbing sets. 

A recently proposed approach |9|, [ 1 1 avoids certain dom- 
inant absorbing sets without compromising code properties 
by carefully selecting the rows/columns of the (SCB) mother 
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matrix. This paper builds upon the cycle consistency matrix 
(CCM) approach ifTTl to analyze SCB codes with column 
weight 5, i.e., codes with five rows of circulant submatrices. 

For an absorbing set to be present in an SCB LDPC 
code, the associated CCM must not be full column-rank. 
Furthermore, analysis of the variable-node graph for a variety 
of absorbing sets reveals that the existence of certain absorbing 
sets is a necessary condition for other (larger) absorbing sets 
to exist. Using these two observations, this paper analyzes the 
smallest absorbing sets in r=5 SCB codes, and systematically 
avoids these absorbing sets by selecting rows from the SCB 
mother matrix that force the associated CCMs to be full 
column-rank. FPGA simulation results confirm that the new 
codes have significantly steeper error-floor slopes . 

As a representative instance of SCB codes, Tanner- 
construction codes [12] with moderate rates are analyzed. It 
is shown that these codes can have improved absorbing set 
spectra by carefully avoiding the smallest absorbing sets with 
suitable code parameter selection. 

Section |ll] describes separable circulant-based (SCB) codes 
and the cycle consistency matrix (CCM). Section III applies 



the CCM approach to analyze three groups of dominant 
absorbing sets in an example family of SCB codes. Section |Hl| 
then selects specific rows from the SCB mother matrix to 
eliminate all the dominant absorbing sets by forcing the 
associated CCMs to be full column-rank. Section [III] also 
analyzes the absorbing set spectrum of Tanner-construction 
codes and provides several good row-selection functions. Sec- 
tion |IV] provides hardware simulation results demonstrating 
the performance improvement obtained with the new codes. 
Section IVl delivers the conclusions. 

II. Definition and Preliminaries 

This section introduces separable, circulant-based (SCB) 
codes and the cycle consistency matrix (CCM) associated with 
absorbing sets in SCB codes. 

A. Circulant-based LDPC codes 

Circulant-based LDPC codes are a family of structured 
regular (r, c) codes where r is the variable-node degree and c 
is the check-node degree. They are constructed as r rows and 
c columns of circulant matrices. They are especially amenable 
to high-throughput hardware implementations |fl3ll . 

The parity-check matrix of circulant-based LDPC codes has 



the following general structure: 
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replacing the difference of j's with the difference of u's , we 
can rewrite the cycle consistency equation as 



where uisapxp circulant matrix. 

A column (row) group is a column (row) of circulant 
matrices. Each variable node has a label (j, k) with j E 
{0, ...,c — 1} being the index of the corresponding column 
group and k E {0, ...,p — 1} identifying the specific column 
within the group. Similarly, each check node has a label 
where i E {0, r — 1} and I E {0, ...,p — 1}. 

Circulant-based LDPC codes include, for example, the con- 
structions in fl2l . Ifl4l and fl5l . The girth can be guaranteed 
to be at least 6 by placing a constraint on the values of the 
submatrix exponent f(i,j) Ell . 

This paper focuses on separable, circulant-based (SCB) 
codes defined as follows: 

Definition 1 (Separable, Circulant -Based (SCB) Code): 
An SCB code is a circulant-based LDPC code with a 
parity-check matrix Hp C f in which f(i,j) is separable, i.e., 
f(i,j)=a(i)-b(j). ' U 

Parity check matrices of SCB codes with the specified 
circulant matrix can be viewed as originating from a common 
SCB mother matrix Hp'j with f m (i, j) = i ■ j. The functions 
a(i) and b(j) effectively specify which rows and columns of 
the mother matrix are selected for the resultant SCB matrix. 
The ranges of a(i) and b(j) are {0, . . . ,p — 1}. 

The SCB structure imposes four conditions [2] on the 
variable and check nodes: (1) bit consistency, (2) check 
consistency, (3) pattern consistency, (4) cycle consistency. 
These conditions are essential to the CCM approach that is 
introduced next. 

B. Absorbing sets and the Cycle Consistency Matrix 

An LDPC code with parity-check matrix H is often viewed 
as a bipartite (Tanner) graph Gh = (V,F,E), where the set V 
represents the variable nodes, the set F represents the check 
nodes, and E corresponds to the edges between variable and 
check nodes. 

For a variable node subset T4 S C V, let G as = (V as , F as , i? as ) 
be the bipartite graph of the edges E RS between the variable 
nodes Va S and their neighboring check nodes F as . Let o(V as ) C 
F- ds be the neighbors of with odd degree (unsatisfied check 
nodes) in G as and e(T4 s ) C F as be the neighbors of V as with 
even degree in G as (satisfied check nodes). An (a, b) absorbing 
set J21 G as = (Va S , F as , E^s) is a Tanner graph with a variable 
nodes, b odd-degree check nodes, and with each variable node 
having strictly fewer odd-degree neighbors than even-degree 
neighbors. 

Suppose there are n variable nodes in the absorbing set. 
Let ji , . . . , j n be the column-group labels of these n nodes in 
the SCB mother matrix. Define u m = j m — ji,m = 2, ...,n 
and u — [u2, ...,u n ]. For each cycle in the absorbing set, by 
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i m )u m = mod p, 



(1) 



where 2t is the cycle length. Note that i m will be different for 
different cycles reflecting the particular cycle trajectories. 

Every cycle in the absorbing set satisfies an equation of the 
form ([TJ. Taken together, these equations produce a matrix 
equation: Mu = mod p, where M ym is the coefficient of 
u m in ([TJ for the yth cycle. 

A key property of M is that Mu = mod p completely 
characterizes the requirement that every cycle in G as satisfies 
([TJ. Even so, it is not necessary for M to include a row for 
every cycle in the absorbing set. 

A cycle need not be included in M if it is a linear 
combination of cycles already included in M. Thus the number 
of rows needed in M is the number of linearly independent 
cycles in G as . Two definitions [16| from graph theory are 
necessary to establish the number of linearly independent 
cycles in G as and hence how many rows are needed for M. 

Definition 2 (Incidence Matrix): For a graph with n ver- 
tices and q edges, the (unoriented) incidence matrix is an n x q 
matrix B with By = 1 if vertex Vi and edge Xj are incident 
and otherwise. ■ 

Definition 3 (Binary Cycle Space): The binary cycle space 
of a graph is the null space of its incidence matrix over GF(2). 
■ 

Any absorbing-set bipartite graph G as can be transformed 
into a graph whose only vertices are V. ds and where two vertices 
are connected iff there is a check node that connects them. We 
call this graph the variable-node (VN) graph of the absorbing 
set. The incidence matrix of the VN graph provides a char- 
acterization of all the cycles in an absorbing set. The number 
of linearly independent cycles in an absorbing set, which is 
the dimension of its binary cycle space, is the size of the null 
space of the incidence matrix B as : Z?bcs = 9 — rank(B as ). 

Having established the number of rows in M, it can be 
formally defined as the Cycle Consistency Matrix: 

Definition 4 (Cycle Consistency Matrix): The cycle consis- 
tency matrix M of an absorbing-set graph G as has |Va S | — 1 
columns and D\, cs rows. The rows of M correspond to 
Dbcs linearly independent cycles in G as . Each row has the 
coefficients of u in ([T| for each of these linearly independent 
cycles. ■ 

Note that M • u = mod p completely characterizes the 
requirement that every cycle in G as satisfies ([TJ. 

The vector u cannot be an all-zero vector because an all- 
zero u indicates that all variable nodes have the same column 
group. This violates the Check Consistency condition, which 
requires that variable nodes sharing a check node have distinct 
column groups. Thus u ^ 0, and a necessary condition for the 
existence of a given absorbing set is that its M does not have 
full column-rank in GF(p). 

If the VN graph of the absorbing set G^ is a sub-graph 
of the VN graph of another absorbing set G^ with the same 



number of variable nodes, then we say the VN graph of the 
absorbing set is extensible. 

Theorem 1: Given a proposed absorbing set graph G as = 
(T4s, F as , E as ), where every variable node is involved in at least 
one cycleQ specified column group labels of the variable nodes 
in Vas in the SCB mother matrix, and specified row-group 
labels of the check nodes in F as in the SCB mother matrix, the 
following are necessary conditions for the proposed absorbing 
set to exist in each daughter SCB LDPC code (with a parity 
check matrix H that includes the specified row and column 
groups of that SCB mother matrix): (1) The CCM for G as does 
not have full column-rank; (2) Variable nodes in V as satisfy the 
Bit Consistency condition and can form a difference vector u 
in the null space of the CCM; and (3) Each check node in 
F as satisfies the Check Consistency condition. Taken together, 
these conditions are also sufficient if the VN graph of this 
absorbing set is not extensible. 



Proof: The proof is provided in Appendix VI-A 



Corollary 1: Given an (ai,bi) absorbing set graph G as = 
(Ks.^as^as) and an (02)62) absorbing set graph G^ = 
(V£,F*,e£), if ai < a 2 and the VN graph of G as is a sub- 
graph of the VN graph of G as , then the existence of G as is a 
necessary condition of the existence of G as . 



Proof: Appendix VI-B provides the proof. 



III. Illustrative case study with r = 5 

The approach in J9l reveals that a careful selection of r 
row-groups from the SCB mother matrix can eliminate certain 
small absorbing sets to improve the error floor. This section 
provides an example with r = 5 (five row groups) that shows 
how to use the new CCM approach to efficiently improve an 
SCB code by analytically avoiding low-weight absorbing sets 
with careful row selections. 

Our example of SCB code design involves three classes of 
SCB codes: (1) Array -based codes ifTTl : the most elementary 
SCB codes in which the first r rows of the SCB mother matrix 
Hp'ftfihj) — i ' j comprise the parity-check matrix. We 
will refer to this class as the elementary array-based (EAB) 
codes. (2) Selected-row (SR) SCB codes: the parity-check 
matrix for these codes is Hp'y,f(i,j) = a(i) ■ j where a(i) 
is called the row-selection function (RSF). (3) Shortened SR 
(SSR) SCB codes: the parity-check matrix for these codes is 
Hpj,f(i,j) — a(i) ■ b(j) where b(j) is called the column- 
selection function (CSF). 

Theorem [T] shows that an absorbing set may be avoided 
either by forcing the associated CCM to be full column rank or 
by precluding u from being in the null space of M. Corollary [T] 
shows that if an absorbing set does not exist, then all absorbing 
sets whose VN graphs contain the VN graph of this absorbing 
set also do not exist. The CCM approach carefully selects the 
RSF and CSF to systematically eliminate small absorbing sets, 
in the order of the size of the VN graph of the absorbing sets. 

Prior results ifTHl prove (4, 8) absorbing sets to be the 
smallest possible for a general r—5 SCB code family and 

' If the variable node degree is at least 2, then each variable node in a 
given absorbing set must be a part of at least one cycle. 





(a) (4, 8) absorbing set 



(b) (5, 9) absorbing set 



Fig. 1: Depiction of absorbing sets. 

to dominate the low BER region iTPTl . Hardware simulation 
results in [9| show that the next two dominant absorbing sets 
in r=5 EAB codes or SR SCB codes are (5,9) and (6,8). 

Careful selection of the RSF via the CCM approach suc- 
cessfully avoids all (4,8), (5,9), and (6,8) absorbing sets. 



Sections [HLA] |III-B| and |III-C| respectively, show that (4,8), 



(5,9), and (6,8) absorbing sets exist in the EAB SCB code 
with r=5 but also that SR SCB codes can systematically 
eliminate these configurations with a carefully selected RSF. 
Section |III-D| provides example RSFs that eliminate small 
absorbing sets with r=5. Section III-E explores the absorbing 
set spectrum of the existing quasi-cyclic LDPC codes with the 
Tanner construction |[T2l . 

A. (4, 8) absorbing sets 

A result in |18| shows that (4,8) absorbing set in Fig. [Ta| is 
the smallest absorbing set in the H p ' P fu j) code family, and |[9] 
shows that certain RSFs eliminate the (4, 8) absorbing sets 
in SR SCB codes. This section uses the CCM approach to 
efficiently establish necessary and sufficient conditions for the 
existence of (4, 8) absorbing sets. 

Lemma 1: det M = mod p is a necessary and sufficient 
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condition for the existence of (4, 8) absorbing sets in Fig 

Proof : Appendix |VI-C| provides the proof. 

Corollaries [2] and [3] restate Lemmas 1 and 2 of 0. The 
CCM approach concisely proves these corollaries in Appendix 
lVLD] 

Corollary 2: (4, 8) absorbing sets exist in EAB codes de- 
scribed by the parity check matrix Hjf „., and their number 
scales as 9(p ). 

Corollary 3: There are no (4, 8) absorbing sets in the SR 
SCB codes described by the parity check matrix H b,p , ^ , for 
prime p large enough with a proper choice of RSF. 

Remark 1: The SR SCB codes avoid (4, 8) absorbing sets 
by carefully choosing RSF such that det M ^ mod p for p 
large enough. One such RSF is [0, 1, 2, 4, 6], and the resulting 
SR SCB codes avoids (4, 8) absorbing sets for prime p > 23. 

B. (5,9) absorbing sets 

Assuming an RSF that avoids the (4, 8) absorbing sets, this 
section proves that the (5, 9) absorbing sets are the smallest 
remaining. The CCM approach shows that the (5, 9) absorbing 
sets always exist in the EAB SCB codes, but are avoided for 



SR SCB codes by some of the RCFs that precluded the (4, 8) 
absorbing sets. 

Lemma 2: Assuming (4, 8) absorbing sets do not exist, 
(5, b) absorbing sets also do not exist for b < 9. 



Proof: Appendix VI-E| provides the proof. ■ 
Corollary 4: Assuming (4, 8) absorbing sets do not exist, 
the (5, 9) absorbing set is the smallest absorbing set. 



Proof: Appendix VI-F provides the proof. ■ 
Lemma 3: det M = mod p is necessary and sufficient 
for the existence of (5, 9) absorbing sets in Fig. [Ib] 

Proof: Appendix |VI-G| provides the proof. ■ 
Appendix VI-H shows that the number of (5, 9) absorbing 



sets in the EAB codes scales as 0(p 3 ). However, with a 
proper choice of RSF, det M/0 mod p for p large enough. 
Therefore we can conclude with the following two corollaries. 

Corollary 5: (5, 9) absorbing sets exist in EAB codes de- 
scribed by the parity check matrix H 'f., and their number 
scales as 0(p 3 ). 

Corollary 6: There are no (5, 9) absorbing sets in the SR 
SCB codes described by the parity check matrix H^,^ ^, for 
prime p large enough with a proper choice of RSF. 

Proof: For SR SCB codes, it is sufficient to select RSF 
such that Q in Appendix VI-G does not evaluate to zero. 
One such example is [0,1,2,4,7], where it is sufficient for 
the prime p > 89 and not be in the set {101, 103, 131, 179}. 
Therefore in the SR SCB code there is no (5, 9) absorbing set 
if p is sufficiently large and an appropriate RCF is chosen. ■ 

C. (6, 8) absorbing sets 

This section considers the (6, 8) absorbing sets, which are 
the smallest remaining after the (4,8) and (5,9) absorbing 
sets. We will investigate the (6, 8) absorbing sets both for 
EAB codes and for SR-SCB codes that preclude the (4,8) 
and (5, 9) absorbing sets. For the six candidate configurations 
of (6, 8) absorbing sets, all satisfied checks have degree 2 
and all unsatisfied checks have degree 1 in the absorbing 
set graph. Combinatorial and consistency arguments show 
that four of these six configurations are not present for p 
sufficiently large in either the EAB code or in SR-SCB 
codes that preclude the (4, 8) and (5, 9) absorbing sets. The 
remaining two configurations have the cardinality 0(p 3 ) in 
the EAB code. However, both of these configurations contain 
a (4, 8) absorbing set as a subset and thus cannot be present 
in SR-SCB codes that preclude the (4, 8) absorbing sets. 
1) (6,8) configuration 1 - check nodes with degree>2: 
Lemma 4: For H p %j and p sufficiently large, there are no 
(6, 8) absorbing sets for which a check node connects to more 
than two variable nodes in the absorbing set graph. 



Proof: Appendix VI-I provides the proof. ■ 
Therefore attention is restricted to the case where all check 
nodes in the absorbing set graph have degree at most 2. 

In a candidate (6, 8) absorbing set, variable nodes can have 
3, 4 or 5 satisfied checks. By the girth constraint, there can 
be at most 2 variable nodes with 5 satisfied checks. Suppose 
there are two such variable nodes. Since there are a total of 8 




Fig. 2: (6, 8) configuration candidate 2. 
unsatisfied checks, the other 4 variable nodes must each have 
3 satisfied and 2 unsatisfied checks. This necessarily implies 
the configuration shown in Fig. [2] which we discuss next. 

2) (6,8) configuration candidate 2 - Fig. [2j Since the 
configuration in Fig. [5] has two overlapping (4, 8) absorbing 
sets, the existence of this (6, 8) absorbing set relies on both 
CCM determinants of these two (4, 8) absorbing sets that share 
an edge. Therefore with a proof similar to that of Lemma [T] 
we can show the following result: 

Lemma 5: The necessary and sufficient condition for the 
existence of (6, 8) absorbing sets in Fig. [2] is det Mi = 
mod p and det M2 = mod p, where Mi and M2 are 
CCMs of the two internal (4, 8) absorbing sets. 

We can also prove that there are p 2 (p — 1) such absorbing 
sets in the EAB codes with any prime p. The proof is shown 



in Appendix VI-J 

Similarly, the following corollaries are consequences of 
Corollary [2] and [3] 

Corollary 7: The (6, 8) absorbing sets in Fig. [51 exist in 
EAB codes described by the parity check matrix H p 'f.j, and 
their number scales as 0(p 3 ). 

Corollary 8: There are no (6, 8) absorbing sets in Fig. [2] 
in the SR SCB codes described by the parity check matrix 
Hp'auyj* f° r prime P large enough with proper choice of RSF. 

Suppose now that there is exactly one variable node in the 
absorbing set having all five checks satisfied. The variable 
nodes in the absorbing set must necessarily be arranged either 
as in Fig. [5] or Fig. [6] 

Similar analysis of (6, 8) configuration candidates 3 to 6 is 
shown in Appendix |VI-K| |VI-L| |VI-M| and |VI-N| respectively. 
D. Non-existence of (4, 8), (5,9) or (6,8) absorbing sets in 
SR SCB codes with a well-chosen RSF 



The following is a consequence of Lemma [T] to Lemma 10 
and Corollary [3] to Corollary 1 1 



Theorem 2: In the EAB SCB code the number of (4,8), 
(5, 9) and (6, 8) absorbing sets scales as Q(p 3 ) whereas in 
an SR SCB code with a well-chosen RSF there are no (4, 8), 
(5, 9) or (6, 8) absorbing sets for sufficiently large p. 

Remark 2: For small p's, (4,8), (5,9) and (6,8) absorbing 
sets cannot be eliminated simultaneously for the SR SCB 
codes. Since there are 3 equivalence classes of SR codes |10|, 
we only need to consider the RSF that contains and 1 . Here 
are a few good RSFs for p > 61 that eliminate all (4, 8), 
(5,9) and (6,8) absorbing sets: (1) [0,1,2,4,17] for p = 67 
(2) [0, 1, 2, 3, 11] for p = 73 (3) [0, 1, 2, 6, 7] for p = 79 (4) 



TABLE I: Error profiles for the EAB SCB (4489, 4158), code 
(top), and the SR SCB code (bottom), n.e. is the number of 
collected errors. 



SNR 


n.e. 


(4,8) 


(5,9) 


(6,8) 


(6,10) 


(7,9) 


(7,11) 


(8,6) 


(8,8) 


(8,10) 


5.6dB 


150 


67 


17 


22 


7 


6 


5 


6 


6 


3 


5.8dB 


139 


83 


18 


16 


6 


5 


1 


3 


3 


1 


6.0dB 


131 


77 


18 


22 


5 


1 


1 


2 


1 


1 


6.2dB 


107 


85 


10 


5 


4 


2 














SNR 


n.e. 


(4,8) 


(5,9) 


(6,8) 


(6,10) 


(7,9) 


(7,11) 


(8,6) 


(8,8) 


(8,10) 


5.6dB 


106 











25 


15 


6 


15 


13 


6 


5.8dB 


140 











35 


29 


14 


16 


6 


8 


6.0dB 


60 











25 


7 


5 


9 


5 


3 



[0,1,2,3,7] for p = 83,97,101,103,107,109,113,127 (5) 
[0, 1, 2, 4, 11] for p = 89 (6) [0, 1, 2, 4, 7] for p > 179. 
E. Absorbing set spectrum in the Tanner construction 

We can easily extend our analysis to the Tanner construction 
in lfT2ll . which is an example of SSR SCB codes. 

Lemma 6: In the Tanner graph corresponding to quasi- 
cyclic LDPC H 5 ' p ~ s in fl/Jl, no (4,8) or (6,8) absorbing 
set exists with parameters selected in Table I of [12|. 



Proof: Proof is shown in Appendix | VI-0 
Remark 3 



This code has a good absorbing set spectrum 
and thus will have low error floor as expected. Moreover, 
the analysis can be easily extended to the quasi-cyclic codes 
constructed in [14], which have the first sub-row and first sub- 
column as identity matrices. 

Note that the codes listed in Table I of lfl2ll are mostly 
moderate-rate codes. However, for higher rates, the Tanner 
construction may introduce smaller absorbing sets ATI . These 
absorbing sets can also be avoided by a carefully chosen CSF 
that precludes the u from the null space of M. 

IV. Results 

In this section we experimentally demonstrate performance 
improvement with well-designed SR codes. In simulations, we 
use 200 iterations and a Q4.2 fixed-point quantization with 4, 
resp. 2, bits to represent integer, resp. fractional, values. We 
simulate sum-product algorithm lfl3l on an FPGA platform. 

We simulated a pair of longer block length codes contrasted 
in Fig. [3] The performance improvement of the SR-SCB 
code is due to the elimination of the (4,8), (5,9) and (6,8) 
absorbing sets with a proper choice of the RSF. (Here the 
SR-SCB code's RSF is [0, 1, 2, 4, 17].) We observe that in the 
error profile as shown in Table [I] the (4,8), (5,9) and (6,8) 
absorbing sets are completely eliminated in the SR SCB code. 

Row selection alone cannot avoid the next smallest ab- 
sorbing sets, which are the (6, 10). Each of these (6, 10) 
absorbing sets exists if and only if the corresponding subset 
(6,4) absorbing set studied in ifTOll exists. As shown in |10|, 
these (6, 4) absorbing sets cannot be eliminated only with row 
selection. However, [ 1 1 also shows that these absorbing sets 
can be precluded by column selection. 

V. Conclusion 
This paper presents a detailed analysis of the absorbing 
set spectrum of a class of LDPC codes based on circulant 
matrices. Using the cycle consistency matrix description of the 



LDPC (p=67, r=5) AWGN fixed point Q4.2 



■ N=4489 EAB code rale=0 9263 
■N=4489 SR code rale=0 9263 



W dB > 

Fig. 3: Performance comparison of the (4489, 4158) EAB SCB 
and SR SCB LDPC codes. 

dominant absorbing sets we characterized code performance 
and provided tools for a systematic code design. Simulation 
results in low FER region support the proposed methodology. 
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VI. Appendix 

A. Proof of Theorem [JJ 

Each of the three conditions has already been shown to be a 
necessary condition for the existence of G as in an SCB. If all 
of these three conditions are satisfied, all the cycles presented 
in the CCM exist in Gh and any linear combination of these 
cycles exists in Gh as well. The only issue is whether the ex- 
isting graphical structures have additional linearly independent 
cycles not required by the CCM. There are only three ways for 
this to happen: (1) a variable node's unsatisfied check node is 
the same as another variable node's unsatisfied check node, or 
(2) a variable node's unsatisfied check node is the same as one 
satisfied check node in the graph, or (3) two of the satisfied 
check nodes are the same. In each of these cases, additional 
edges extend the VN graph. For this to be possible, the original 
VN graph must be extensible as defined above. Thus if the VN 
graph is not extensible, the above constructed solution fully 
describes the existence of the proposed absorbing set. This 
concludes that G as is present in Gh- 

B. Proof of Corollary [7J 

Suppose the CCMs of G\ & and G% s are Mi and M 2 
respectively. If the VN graph of G as is a sub-graph of the 
VN graph of G as , the independent cycles of G as will also be 
independent in G^ and thus Mi could be a sub-matrix of M 2 : 



M 2 = 



Mi 
A 




B 



(2) 



Therefore if there exists a valid u 2 such that M2U2 = 
mod p, the first a x elements would also be a valid u x such 
that M1U1 = mod p. This concludes that the existence of 
G as is a necessary condition of the existence of G as . 

C. Proof of Lemma [JJ 

Since the VN graph of the (4, 8) absorbing sets is a fully 
connected graph, it is not extensible without introducing a 
length-4 cycle in the corresponding bipartite graph. According 
to Theorem [T] det M = mod p is a necessary condition 
since it implies the M is not full column-rank, and it is a 
sufficient condition if bit consistency and check consistency 
are both satisfied. We can verify both of these conditions by 
carefully constructing a solution in the null space of M. 
Since the binary cycle space for Fig. la has dimension 3, 



we construct the following CCM by selecting the following 3 
linearly independent cycles: v\ — v 2 — V3, v\ — v 2 — V4, v\ — 

v 3 - v 4 : 

ii — 12 12 — is 
ii — ia 16 — «4 
is — 13 «3 — ii 



M : 



(3) 



If det M = mod p, there exists a non-zero solution to 
M • u = mod p, where u = [u 2 , 1*3,^4] T . Without 
loss of generality, suppose u 2 7^ 0. With check consistency, 

i\ — 12 7^ 0,«2 — H 7^ 0, and i\—i 2 7^ ii — H- Thus 113 7^ and 



u 2 7^ u 3- Similarly u 4 7^ 0, u 2 7^ U4, and u 3 7^ U4. Then, for a 
fixed ji, we can find j 2 , jz, and j'4 without contradiction to the 
bit consistency. With any specific k\ values, we can find (4, 8) 
absorbing sets in the code. Therefore det M = mod p is a 
sufficient condition for the existence of (4, 8) absorbing sets. 

D. Proof of Corollary [2] and Corollary [i] 

From the proof in Appendix VI-C| det M = mod p is a 
necessary and sufficient condition for the existence of (4, 8) 
absorbing sets. 

Under the bit consistency and check consistency, there 
are only two possible non-isomorphic check labelings 
0: (ii,i 2 ,«3,*4,«5>«6) to be either (x,y,x,y, z,w) for 
{x, y, z, w} C {0, 1, 2, 3, 4} (assignment 1) or (x, t, w, y, z, z) 
(assignment 2). Thus det M = mod p implies the follow- 
ing necessary and sufficient conditions for the existence of a 
(4, 8) absorbing set result: 

(z — x)(w — y) + (z — y)(w — x) = mod p , (4) 

for assignment 1, and 

(z— w)(x— t)(y— z) — {y— w)(x—z)(z— t) = mod p , (5) 

for assignment 2. 

For EAB codes {x,y,z,w,t} — {0,1,2,3,4}, in 
the former case, in fact one can show there are no 
solution sets for prime p large enough (p > 17). For 
the latter case there are 8 solution sets (x,y,z,w,t) E 
{(4, 3, 2, 0, 1), (4, 1, 2, 0, 3), (3, 4, 2, 1, 0), (3, 0, 2, 1, 4), 
(1, 4, 2, 3, 0), (1, 0, 2, 3, 4), (0, 3, 2, 4, 1), (0, 1, 2, 4, 3)} that 
always evaluate to zero on the left-hand side of equation (PJ. 
These numerical solutions are in fact symmetric so that once 
the labels of the check nodes are selected (cf. Fig. [la) , the 
variable node labels (pairs through (^4,^4)) can be 

selected in 8(p 3 ) ways, thereby completely characterizing 
the absorbing set of interest. 

We say that an absorbing set is fully absorbing set if, in 
addition, all variable nodes outside the absorbing set have 
more satisfied than unsatisfied checks. 

Then the above (4, 8) absorbing set is always a (4, 8) fully 
absorbing set since otherwise there would exist a variable node 
(j5j^5) outside the absorbing set incident to at least three 
of the checks labeled 17 through 114. Such a configuration 
would either violate the girth constraint, or it would imply 
the existence of a new configuration spanning four variable 
nodes (these being the node k^) and three variable nodes 
from the starting (4, 8) absorbing set). These four variable 
nodes would necessarily be connected such that their common 
constraint is given in Q, previously shown to not hold for 
large enough p. 

To avoid this absorbing set we need to force the determinant 
of the CCM to be nonzero. It is sufficient to assign values 
to a(i) such that the selected labels for the check nodes 
i\ through i§ do not satisfy equations Q and {3). One 
such example is (i,a(i)) e {(0, 0), (1, 1), (2, 2), (3, 4), (4, 6)}. 
Therefore in the SR SCB code there is no (4, 8) absorbing set 
for large enough p. 



E. Proof of Lemma [2] 

Since the total edge number is odd and the number of 
edges that go to satisfied check nodes is even, b can only 
be odd. Thus if b < 9, b can only choose values from 
{1,3,5,7}. For these b values, the corresponding VN graph 
will always contain the VN graph of (4, 8) absorbing sets 
without introducing a length-4 cycle in the bipartite graph. 
By Corollary [I] the (5, b) absorbing sets with b < 9 do not 
exist if (4, 8) absorbing sets are absent. 

F. Proof of Corollary [5] 

Since the number of edges that go to unsatisfied check nodes 
is at most 10 and b is odd, the only possible (5, b) absorbing 
set is (5,9) absorbing set in the absence of (4,8) absorbing 
sets. Since the only possible configuration of a (5, 9) absorbing 
set is one variable node that has 4 satisfied check nodes and 



lb 



four variable nodes that have 3 satisfied check nodes. Fig 
depicts this configuration, which has a VN graph that does 
not contain the VN graph of (4, 8) as a sub-graph. Therefore 
we can use Lemma [3] to build the necessary and sufficient 
condition of this absorbing set. 

G. Proof of Lemma [3] 



The binary cycle space for Fig. lb has dimension 4. We 
construct the following CCM by selecting the following lin- 
early independent cycles: v\ — v 2 — v%, v\ — v 2 — v$, Vi — V3 — 

v 4 ,Vi -V4- v 5 : 



M 



Zl — *5 «5 — *2 

ii — is 

12 — is 






is — 14 
i6 — i'i 

is — i7 %t — ii 



(6) 



Similarly to the proof of Lemma [T] we can show that 
det M = mod p implies the following necessary and 
sufficient conditions for the existence of a (5,9) absorbing 
set, where 



(7) 



detM - i 5 )(is - u)(l2 - ia)(i'i - ii) 

- (ii - is)(is - «2)(«6 - «3)(«7 - U) mod p 

H. Number of (5, 9) absorbing sets in EAB codes 

Under the bit consistency and check consistency, 
there are 5 possible non-isomorphic check labelings: 

*2, *3, *4, *5 5 *6, *s) to be (x,y, z,w, z,w,x,y), 
(x, y, z, w, t, w, x, y), (x, y, z, w, t, w, x, z), 

(x,y, z,w,t, x,t,y), or (x,y, z,w,t,w,t.y). The EAB 
codes have {x,y,z,w,t} C {0,1,2,3,4}. For the 
4th case, there are 8 solution sets (x,y,z,w,t) E 
{(4, 0, 1, 3, 2), (4, 0, 3, 1, 2), (3, 1, 4, 0, 2), (3, 1, 0, 4, 2), 
(1, 3, 4, 0, 2), (1, 3, 0, 4, 2), (0, 4, 1, 3, 2), (0, 4, 3, 1, 2)} that 
always evaluate to zero on the left-hand side of equation f7i. 
Once the labels of the check nodes are selected (cf. Fig. [lb I, 
the variable node labels (pairs through (^5,^5)) can 

be selected in 0(p 3 ) ways, thereby completely characterizing 
the absorbing set of interest. 




Fig. 4: (6,8) configuration candidate 1. 



/. Proof of Lemma |4] 

For those (6, 8) absorbing sets with the property that a check 
node connects to more than 2 variable nodes in the absorbing 
set, the corresponding VN graph will always contain the VN 
graph of (4, 8) absorbing sets without introducing a length-4 
cycle in the bipartite graph. One such example is shown in 
Fig. [4] Since this structure contains several (4, 8) absorbing 
sets, we establish a CCM for each (4, 8) absorbing set. If 
the configuration in Fig. |4] were to exist, the CCMs of all 
(4, 8) absorbing sets within this configuration would have zero 
determinant, resulting in a set of conditions which cannot be 
satisfied simultaneously for p large enough. 

J. Proof of Lemma [5] 

We may start with the substructure spanning 4 variable 
nodes (j x ,k\), (j 2 ,k 2 ), {j 5 ,h) and (j e ,h). Recall that out 
of six check nodes shared by these variable nodes, exactly 
two have the same label (it being label "2"), and that there 
is p 2 (p — 1) ways of assigning values of the variable nodes 
and the check nodes in this substructure. By symmetry of 
the configuration it suffices to consider the case when this 
repeated label is the one corresponding to 15 and iio and when 
this repeated label corresponds to some other pair of parallel 
edges. The latter case is not possible since the variable nodes 
(ii)fci), (j2,&2), (.73,^3) and (j4,fc4) themselves constitute 
a (4, 8) absorbing set and there the propagation of the check 
labels through the proposed configuration necessarily violates 
the check label constraints. In the former case, by the bit 
consistency constraints, the labeling of the checks incident 
to (^3,^:3) and (j4,A;4) is unique for each of the nodes 
(without assuming these two nodes themselves share an edge). 
Moreover, for the independently selected values of (j'3, k^) 
and (^'4,^4), we show by the pattern consistency constraint 
that they indeed have a common check iio, itself labeled "2". 
As a result, once the values (ji, fei), (^2,^2), (^5,^5) and 
(j6j ^6) and their shared checks are pinned down - which can 
be done in p 2 (p — 1) ways - the rest of the proposed (6, 8) 
configuration follows uniquely. As a result, the cardinality of 
(6, 8) absorbing sets of the type in Fig. [2] is p 2 (p — 1). 



Fig. 5: (6, 8) configuration candidate 3. 



Fig. 6: (6, 8) configuration candidate 4. 
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K. (6, 8) configuration candidate 3 - Fig. [5] 

Lemma 7: In the Tanner graph corresponding to . 
for the EAB SCB and for the SR SCB code there are no (6, 8 
absorbing sets for p large enough of the type shown in Fig. [5] 

Proof: Without loss of generality we may assign check node 
labels for the checks emanating from the variable node k^) 
as follows: h = x,t2 = y> *3 = z, *4 = w, and 25 = t, where 
x, y, z, w, t are the five distinct check labels. 

The binary cycle space for Fig.|5]has dimension 6. We con- 
struct the following CCM by selecting the following linearly 
independent cycles: v% — v% — V3, v\ — U3 — 1>6, v\ — V2 — «5, v\ — 

v 3 - V 4 ,Vl 



M 



(8) 

The rank of the matrix is at most 5 so if fact we may consider 
the top-left 5 by 5 submatrix (call it B). If the matrix B is full 
rank, the only solution is j\ = 22 = 33 = 3a = 35 = 36- Hence 
det(i3) = is necessary for the existence of the absorbing sets 
of this type. Such condition can be expressed as 

- (in - z)[-(x - i e )(i 9 - t)(i 7 - w)(y - ho) 

+ (x- i 7 ){i 6 - t)(y - ig)(ho - w)} = mod p. 

Also consider the bottom-right 4 by 4 submatrix (call it A). 
If the matrix A is full rank, the whole matrix will be rank 5, 
which only have one solution as ji = 32 = 33 = 3a = is = 
3q. Hence det(A) = is necessary for the existence of the 
absorbing sets of this type. Such condition can be expressed 
as 

-(x - i 7 ){i w - w)(y - i n )(i 8 - z) 
+ (x - i 8 )(i 7 - w)(i n - z)(y - ho) = mod p. 

For the values of ie>i7>i9,ho an d in in the set described 
by f(i,j) for both the EAB SCB and for the SR SCB codes, 
and such that the labels meeting at the same vertex are distinct 
(see Pig. [5j, the equation |9| and |9]) evaluate to zero for only 



(9) 



(10) 



finite number of values of the parameter p. For a(i) = i (EAB 
SCB code) det(i?) 7^ for p > 23. For the carefully designed 
SR SCB code, the equation |9]l and (|9]l also evaluate to zero 
for only finite number of values of the parameter p. Thus for 
the EAB SCB code and for the SR SCB code there are no 
(6, 8) absorbing sets for p large enough of the type shown in 
Fig.0 ^ " ■ 

Corollary 9: In the Tanner graph corresponding to 
(^'^) a bsorbing sets of the type shown in Fig. |5j 
exist if and only if det(A) = and det(B) = 0. 

Proof: Since the rank of a matrix is lower-bounded by the 
rank of its submatrix, and M has a submatrix 



i 9 -t 

x — ij h — w 

hi-z 

iio-«; 



(11) 



which is a full-rank 4 by 4 matrix, the rank of M is no less 
than 4. 

If absorbing sets of the type in Fig. [5] exist, M is not full 
column-rank, otherwise the equation only has one solution as 
ji = h = h = k = 35 = 36- With the fact that rank(M) > 
4 and rank(M) < 5, we can have rank(M) = 4, which 
implies det(A) = 0. (Otherwise M will have rank 5 with the 
additional column.) Without loss of generality, we can consider 
the bottom row of A is dependent of the other 3 rows, and 
thus this row is redundant. Then rank(_B) — rank(M) = 4, 
which implies det(i?) = 0. 

If det(A) — and det(£?) = 0, without loss of generality, 
we consider the bottom row of A is dependent of the other 
3 rows and thus rank(S) = rank(M). Due to det(B) = 0, 
M is not full rank and rank(M) < 5. With the fact that 
rank(M) > 4, we have rank(M) = 4 and there exist a non- 
zero u such that equation ( p~2] > is satisfied. Each element of u 
is non-zero, otherwise the u will be a zero vector. 

This completes the proof of the corollary. ■ 

L. (6, 8) configuration candidate 4 - Fig. [6] 

Lemma 8: In the EAB codes corresponding to H p '?p there 
are no (6, 8) absorbing sets with the topology shown in Fig. [6] 
for p large enough. 




belled by either (i x , i 2 , k, k, k, k) 

{k,k,k,k,k,k) = (y,x,z,w,t,z). 



(y,z,x,z,t,w) or 



h2 iu 

Fig. 7: (6, 8) configuration candidate 5. 

Proof: The binary cycle space for Fig. [6] has dimension 6. 
We construct the following CCM by selecting the following 
linearly independent cycles: v% — v 2 — v 3 , V\ — v 2 — V4, v% — 

w 3 - v 4l vi -v 3 - v 6 ,vi - « 4 - u 5 ,wi - «5 - w 6 : 
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assignment 



It suffices to 

k,k,k,U,k,k 
(k,k,k,h,h,k) 

(k,k,k,U,k,k) 

Analogously to proof of Lemma [7j, the M is not full 
column-rank if det 4/0 mod p and det B ^ mod p, 
where A is the left-top 3 by 3 submatrix of M and B is the 
right-bottom 4 by 4 submatrix of M. This constraint cannot 
be satisfied for p > 41 for i\ to in taking values in the set 
{0, 1, 2, 3,4} and such that the bit consistency constraints are 
satisfied for both labellings. ■ 

Note that the variable-node of the configuration shown in 
Fig. [6] is a subgraph of the VN graph of (4, 8) absorbing sets. 
The following again is an easy consequence of Corollary [3] 

Corollary 10: If (4, 8) absorbing sets are absent, the con- 
figuration shown in Fig.|6]is not possible in the SR SCB code. 

M. (6, 8) configuration candidate 5 - Fig. [7] 

In the remainder we consider the case when no variable 
nodes in the absorbing set has all five satisfied checks. This 
constraint implies configurations shown in Fig. [7] and Fig. [8] 

Lemma 9: In the EAB codes corresponding to H^f^ there 
are 0(p 3 ) (6, 8) absorbing sets of the type shown in Fig. [7] 
for p large enough. 

Proof: Note that by the property of (4, 8) absorbing 
sets satisfied check nodes in the substructure spanning bit 
nodes (ji,fci), {h,k 2 ), (j 3 ,k 3 ), and {j4,k 4 ) can be la- 



Using the technique of Section [Tl-B| we construct the CCM 
for this configuration. The binary cycle space for Fig. [7] has 
dimension 5. We construct the following CCM by selecting 
the following linearly independent cycles: v\ — v 2 — v 3 , v\ — 

v 2 ~ U4, vi - v 2 - v 5 , vi - v 3 - v e - t>4, V1 — V5 — V6 — v 4 
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The determinant of the CCM is 



(13) 



detM 

=[— (z — w)(i 5 — i4)((w — iio)(iio — in) — (k — «io)(iio — «<>)) 
+ (i 9 — w)(w — y)(i 5 — ii)(iio — in)](x — t) 
— (x - i 5 )(t - z)(w - y)(iio - i9)(in - £4) mod p 

(14) 

In fact this determinant always evaluates to zero for every 
value of p for 6 non-isomorphic edge labellings. For each such 
labelling, once say the values of ji , ki and j 2 are selected 
(which can be done in p 2 (p— 1) ways), the rest of values in the 
configuration follows uniquely. Therefore there are 6p 2 (p— 1) 
such absorbing sets. 

For certain small values of p, p € {2,3,7, 11, 13,31,47}, 
( fl"4| ) has additional solutions, raising the total number of 
solutions to 8p 2 (p — 1). Nonetheless, the scaling <d(p 3 ) of 
the cardinality of the absorbing sets still holds. 

In the latter case we likewise establish a matrix relating the 
labels of the check nodes (i 7 , i s , k,iio, in in Fig. [6]) incident 
to bit nodes (j'5, k$) and (jg, k$) to adjacent checks. Again, by 
imposing the bit consistency conditions, we conclude that the 
above constraint cannot hold for p > 23 and therefore such a 
labelled configuration is in fact not possible for large enough 
p. U 

As in the previous case, the variable-node of the configura- 
tion shown in Fig. [7] is a subgraph of the VN graph of (4, 8) 
absorbing sets. The following again is an easy consequence of 
Corollary [3] 

Corollary 11: If (4, 8) absorbing sets are absent, the con- 
figuration shown in Fig.|7]is not possible in the SR SCB code. 

N. (6, 8) configuration candidate 6 - Fig. [s] 

The last configuration we consider is the one in Fig. [8] 
Lemma 10: In the Tanner graph corresponding to H p f^ ^ 
there are no (6, 8) absorbing sets for p large enough of the 
type shown in Fig. [8] in neither the EAB SCB nor in the SR 
SCB code. 

Proof: The binary cycle space for Fig. [8] has dimension 
6. We construct the CCM by selecting the following linearly 
independent cycles : v\ — v 2 — v 3l vi— v 2 — 174 , ui — v 2 — V5,v± — 




Fig. 8: (6, 8) configuration candidate 6. 
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(15) 



The inspection of M reveals that rank(M) = 4 if and only if 

det(A) = and det(B) = 0, where 
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Their determinants are 

det(A) = (i 8 - i 2 )[— (ii - i4)(ia - is)(i4 - is) 

+ (i 4 — iio)((i3 — ie)(i3 — U) - («2 — i3)(«5 — 13))] 

(18) 

det(B) = (ii — iio)[— (is — *2)(ig — hi)(ii — i3)(«6 — 17) 
+(«2 — ii)((ie — ii)(ie — »7)(i9 — in) 
-{ia - ie)((i9 - ie)(i9 - ill) - (h - M)){is - ig)))] 

(19) 



Then in the similar manner of the proof of lemma [7] we 
can show that the determinant of the corresponding matrix 
evaluates to zero only in finitely many choices for p for either 
selection of a(i). In particular for a(i) = i it suffices for p to 
be > 29 and p 7^ 41 for the configuration not to exist. ■ 

O. Proof of Lemma [6] 

For the codes with girth greater than 6, both (4, 8) and (6, 8) 
absorbing sets do not exist since these two sets contain cycle- 
6. For the codes with girth = 6, we take p = 31, c = 5, r = 6 
as an example, with b = 2, a — 6. Thus the H p 'f(ij) is a 
sub-matrix of array code -^'/(j j) w ith f(hj) — Tn(i) ■ j 
for (i,m(i)) e {(0,1), (1,2)^(2, 4), (3, 8), (4, 16)}. We 
set up the system of equations as before. For (4, 8) 

absorbing sets, the only possible labelling for p = 31 is 

(ii, *2> *3j Hi *5 , k) = (x,t,w,y, z, z) and there are five 
non-isomorphic solutions to the equation (z ~ w)(x — t)(y — 
z ) ~ (y — w)(x — z)(z — t) = mod p: (x,y, z,w,t) = 
(2, 4, 1, 16, 8), (1, 8, 2, 4, 16), (1, 2, 4, 16, 8), (1, 4, 8, 2, 16), 
(1, 2, 16, 8, 4). Each of the solutions corresponds to the matrix 



R 



z — x x — t t — z 
y — x x — z z — y 



(20) 



such that R(ji,j2,j3,ji) T = mod p. Suppose the null 
space of each matrix is A^,l < i < 5. For any (4,8) 
absorbing set, (ji, J2, 73, J4) should be in U iVj. Denote 

1<2<5 

Y = {1,6,5,30,25,26}, which is the index of the sub- 
columns in the quasi-cyclic code, (ji, 72, J3, J4) should also 
be in u (Y(i),Y(j),Y(k),Y(l)). However, in this case, 

{ u fflf),Y(j),Y(k),Y(l))} n { U Ni} = 

(4, 8) absorbing sets do not exist, and consequently 
absorbing sets in Figs. [2] [6] and |7] also do not exist. 

Similarly, the null space of equation ( fT2| and 
in lemma [10] and [7] does not have intersection 

U (Y(i),Y(j),Y(k),Y(l)), which eliminates the possibil- 

ity of existence of (6, 8) absorbing sets in Figj5] and [8] 



Thus 



G3 

with 



